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We investigate the ionic Hubbard model on a triangular lattice at three-quarters filling. This 
model displays a subtle interplay between metallic and insulating phases and between charge and 
magnetic order. We find crossovers between Mott, charge transfer and covalent insulators and 
magnetic order with large moments that persist even when the charge transfer is weak. We discuss 
our findings in the context of recent experiments on the layered cobaltates A0.5C0O2 (A=K, Na). 
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The competition between metallic and insulating 
states in strongly correlated materials leads to many 
novel behaviours. The Mott insulator occurs when a sin- 
gle band is half-filled and the on-site Coulomb repulsion, 
U , is much larger than hopping integral, t. A menagerie 
of strongly correlated states is found when a system is 
driven away from the Mott insulating state, either by 
doping, as in the cuprates [l[ , or reducing U /t, as in the 
organics Q ■ Geometric frustration causes yet more novel 
physics in Mott systems 0] ■ Therefore the observation of 
strongly correlated phases in the triangular lattice com- 
pounds A0.5C0O2, where A is K or Na [3J, has created 
intense interest. 

An important model for investigating insulating states 
in correlated materials is the ionic Hubbard model. On 
a half-filled square lattice this model displays a crossover 
between Mott and band insulating states which has been 
analyzed with quantum Monte Carlo (QMC) Q, dynam- 
ical mean field theory (DMFT) and its cluster extensions 
Q. However, except for the case of one dimension 
this model has not been studied away from half-filling 
and/or on geometrically frustrated lattices. 

In this Letter we study the ionic Hubbard model on 
a triangular lattice at three-quarter filling. This Hamil- 
tonian displays a subtle interplay between metallic and 
insulating phases and charge and magnetic order. It has 
regimes analogous to Mott, charge transfer @, and cova- 
lent insulators . The study of this model is motivated 
in part by our recent proposal (Io| that it is an effective 
low-energy Hamiltonian for Na^OoO^, at values of x at 
which ordering of the sodium ions occurs. 

The Hamiltonian for the ionic Hubbard model is 



(i) 



where c^J anihilates (creates an electron with spin a at 
site i, t is the hopping integral, U is the effective Coulomb 
repulsion between electrons on the same site, and ej is 



a the site energy. We specialise to the case with two 
sublattices, A (ej = A/2) and B (e< = —A/2), consisting 
of alternating rows, with different site energies on the two 
sublattices (c.f., Fig. 15 of Ref. [k]). This is the lattice 
relevant to A0.5C0O2 where the difference in site energies 
results from the ordering of the A-atoms [ll|, [IH, [l3[ • 

Two limits of model ([T]) at 3/4- filling may be easily un- 
derstood. For non-interacting electrons, U — 0, a metal- 
lic state occurs for all A as at least one-band crosses the 
Fermi energy. In the atomic limit t = 0, and U > A 
one expects a charge transfer insulator with a charge gap 
of about A whereas for U < A a Mott insulator with 
charge gap of U occurs. However, realistic parametriza- 
tion of AzCoO^ materials imply U ^> A and A ~ \t\ 
[Til ]: we will show below that in this parameter regime 
the model show very different behaviour from either of 
the limits discussed above. This interesting regime needs 
to be analyzed using non-perturbative and/or numeri- 
cal techniques. Thus, we have performed Lanczos diag- 
onalization calculations on 18 site clusters with periodic 
boundary conditions. 

In Fig. [1] we plot the charge transfer, ns — tia as a 
function of A/|i| for several values of U. We also plot 
iib — tia in two analytically tractable limits: the non- 
interacting limit, U = [151 ]; and the strong coupling 
limit U A ^> \t\ [l6j]. Several interesting effects can be 
observed in this calculation. Firstly, the sign of t strongly 
effects the degree of charge transfer on the triangular 
lattice. Secondly, charge transfer depends only weakly on 
U. Thirdly, regardless of the sign of t or the magnitude 
of U, the charge transfer increases rather slowly as A 
increases. 

The charge gap, i.e., the difference in the chemical po- 
tentials for electrons and holes, is A c = Eq(N + 1) + 
E (N - 1) - 2E (N), where E (N) is the ground state 
energy for N electrons. We plot the variation of A c with 
A for various values of U in Fig. [2] A c vanishes for 
[7 = 0, however finite size effects mean that we cannot 
accurately calculate A c for small A. A c = A for t = 
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FIG. 1: (Color online) Variation of the charge transfer with 
A for different Coulomb repulsion energies. Lanczos results 
on an 18-site cluster for: U/\t\ — (stars), 15 (open symbols) 
and 100 (filled symbols) are shown. The U = results are 
in excellent agreement with the infinite lattice tight-binding 
(tb) result (dashed-dotted line) [l5||. Strong coupling results 
(dashed lines) are also shown for comparison. 
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FIG. 2: (Color online) Variation of the charge gap, A c , with 
A for U/\t\ — 15 and 100 (the charge gap is zero for U = 0). 
The tendency towards an insulating state is greater for t < 
than for t > 0. Dashed lines show the strong coupling limit: 
U >> A >> t. Note that although strong correlations are 
essential for creation of the charge gap they are not required 
for the charge transfer, c.f.. Fig. [1] Note that, for U >> A 
with U large, the charge gap is robust against the value of U. 



and U ^S> A; this result is reminiscent of a charge trans- 
fer insulator Q. Both perturbative [13] and numerical 
results show that the charge gap depends on the sign of 
t due to the different magnetic and electronic properties 
arising from the geometrical frustration of the triangu- 
lar lattice. In contrast, on a square lattice, A c , does not 
depend on the sign of t. 

In the limit, A ^> U 3> \t\, the A and B sublattices are 
well separated in energy; the B sites are doubly occupied 
(i.e., the B-sublattice is a band insulator) and the A sub- 



lattice is half-filled and hence becomes a Mott insulator. 
If there were no hybridisation between that chains, one 
would find a metallic state for any finite charge trans- 
fer from the B-sites to the A-sites (self doping), even for 
U ^> \t\ as the A-chains are now electron-doped Mott 
insulators and the B-chains are hole-doped band insula- 
tors. However, Fig. [2] shows that the insulating regime 
of the model extends far beyond the well understood 
iib — riA = 1 regime. This is because the real space in- 
terpretation is incorrect as hybridization between A and 
B chains is substantial. For \t\ ~ A <C U the system can 
remain insulating with a small gap [0(i)]. This state is 
analogous to a covalent insulator [9( . 

One expects that for A = the ground state is metal- 
lic as there the system is 3/4- filled. However, a small but 
finite A = + leads to a strongly nested Fermi surface for 
t > whereas for t < the Fermi surface rather feature- 
less. Thus, rather different behaviors might be expected 
for different signs of t even at weak coupling. At large 
U our exact diagonalization results suggest that a gap 
may be present even for a small value of A/t. However 
finite-size effects, inherent to the method, mean that it 
is not possible to resolve whether a gap opens at A = 
or at some finite value of A. 

To test this covalent insulator interpretation in the 
A ~ \t\ and large U regime we have also calculated the 
spectral density, A(ui), c.f., Fig. [3l There are three dis- 
tinct contributions to the A(uj): at low energies there is 
a lower Hubbard band; just below the chemical poten- 
tial (u> = fx) is a weakly correlated band; and just above 
uj = fi is the upper Hubbard band. Furthermore, the 
large energy separation, much larger than the expected 
U = 15|f|, between the lower and upper Hubbard bands 
is due to an upward (downward) shift of the upper (lower) 
Hubbard bands due to the strong hybridization. In con- 
trast, in the strong coupling limit A(u>) has a much larger 
gap, 0(A), between the contributions from the weakly 
correlated band and the upper Hubbard band. 

The magnetic moment associated with the possible an- 
tiferromagnetism, m v = (3{Sf S^}) 1 / 2 , where v = A or 
B and S? = ^(n^ — riijj, is evaluated between two next- 
nearest neighbors on the v sublattice at the center of clus- 
ter (to reduce finite size effects [Hj])- Fig. [4] shows that 
rriA increases with A and is substantially enhanced by U, 
whereas tub is always small. This is in marked contrast 
to a spin density wave, as predicted by Hartree-Fock cal- 
culations where the magnetic moment is far smaller than 
that experimentally observed (lfjj . 

We now turn to discuss the consequences of our re- 
sults for understanding experiments; for simplicity and 
concreteness we focus on Na^CoC^. The x = 0.5 mate- 
rials have remarkably different properties from those on 



other values of x [20l. l2l|. Above 51 K the intralayer re- 
sistivity of Nao.5Co02 is weakly temperature dependent 
with values of a few mflcm [20( characteristic of a bad 
metal [lQ]. Below 51 K the resistivity increases, consis- 
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FIG. 3: (Color online) Evolution from a charge transfer insu- 
lator to a covalent insulator. The energy dependence of the 
spectral density, A(lj), is shown for two different parameter 
regimes of the model with t < and U = 15[t[. The spectral 
density in the upper panel (A = 10|t|) is that characteris- 
tic of a charge transfer insulator [j| there is a single weakly 
correlated band largely associated with the B sites and ly- 
ing between lower (LHB) and upper (UHB) Hubbard bands 
(separated by ~ U) that are largely associated with the A 
sites and A c ~ A. In contrast, the lower panel (A = \t\) 
shows a spectral density characteristic of a covalent insulator 
(jj: there are only small difference between A and B sites, the 
separation of the LHB and UHB is > U, and A c ~ \t\. 




FIG. 4: (Color online) The magnetic moment as a function 
of A for t < and various U. The magnetic moment on the 
A sites, tua, (main panel) is strongly enhanced when U 3> 
\t\. The inset shows the moment on the B sublattice. tub 
is much smaller, reduced by A, and only weakly dependent 
on U. These results demonstrate that, for large U and A •> 
|t|, the electrons on the A-sublattice are much more strongly 
correlated than those on the B-sublattice despite the small 
charge transfer (see Fig. [TJ. 



tent with a small gap opening (~10 meV) [20j. Thus 
a (bad) metal-insulator transition occurs at 51 K. The 
insulating state of Nao.sCoC^ has a number of counter- 
intuitive properties, not the least of which is the absence 
of strong charge ordering. NMR observes no c harg e or- 
dering up to a resolution of tlb — nA < 0-4 
while neutron crystallography suggest ii B — n,\ 



m 123, 

~ 0.12 

ll| . Thus the insulating state is not the simple charge- 



transfer-like state predicted by (P) in the strong cou- 
pling limit. Nao.5Co02 develops a commensurate mag- 
netic order below 88 K [22|, Hj] . A large magnetic mo- 
ment [m — 0.26(2)^ per magnetic Co ion] is observed 
in spite of the weak charge order [note that classically 
m < (riB — nA)fJ-B/2}. Above 100 K the optical conduc- 
tivity [2l| shows no evidence of a Drude peak, consistent 
with a bad metal. In the insulating phase spectral weight 
is lost below ^10 meV, consistent with the gap seen in the 
dc conductivity and a peak emerges at ~20 meV, which is 
too sharp and too low energy to correspond to a Hubbard 
band [2l|. ARPES shows that the highest ener gy o ccu- 
pied states are ~10 meV below the Fermi energy 25j. No 
equivalent insulating state is seen in the misfit cobaltates 
(221 ] . which supports the contention that Na-ordering is 
vital for understanding the insulating state. 

Various theories have been proposed to explain these 
intriguing experiments. Lee et al. j2o| have performed 
LDA+U calculations, which include Na-ordering, but 
not strong correlations. Other groups [2jj have studied 
strongly correlated models that include the Coulomb in- 
teraction with neighbouring sites, but neglect the effects 
of Na-ordering. Marianetti and Kotli ar 1 281] have also 
studied the Hamiltonians proposed in [101 ] for x = 0.3 
and 0.7. 

In order to compare our results with experiments on 
Nao.5Co02 we need to restrict ourselves to the relevant 
parameter values: t < and \t\ ~ A <C U 14]. This cor- 
responds with the regime of the three quarters filled ionic 
Hubbard model that is both the most interesting and the 
most difficult to study via exact diagonalisation because 
of the deleterious finite size effects. Nevertheless we pro- 
pose that in Nao.sCo02 the insulating state is analogous 
to a covalent insulator. This explains a wide range of ex- 
periments. The peak observed at u> ~ 30 meV in the op- 
tical conductivity [2l[ , is interpreted as the transfer of an 
electron from the weakly correlated band to form a dou- 
blon in the strongly correlated band. The weak charge 
transfer (ng — ua = 0.1 — 0.3; c.f., Fig. [T|) is caused by 
the strong hybridisation between the A and B sublattices 
and is consistent with the value extracted from crystal- 
lographic experiments (0.12 



llj]) and the bounds from 



NMR (<0.4 (2J). The large moment (0.1-0. 2^ B ; c.f., Fig. 
[4]) is comparable to the moment found by neutron scat- 
tering (0.26/^s [24]) and results from the electrons in the 
strongly correlated band, i.e., the single spin hybridised 
between the A and B sublattices. Finite size effects mean 
that we cannot accurately calculate the charge gap in this 
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regime. However, we propose that the experimental sys- 
tem corresponds to a parameter range where the gap is 
small, A r < 0(|i|), consistent with the gap, ~7-10 meV 
[20I [2^ , seen experimentally in ARPES and resistivity. 
This is consistent with the expectation that A c — > as 
A/|i| — ► 0. Accurately calculating A c for small A/|t| 
and large U, and hence further testing our hypothesis, 
therefore remains an important theoretical challenge. 
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